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1. Introduction 

Matrix models |jl|] are an efficient way to describe random surfaces, noncritical string 
tfieory, and perliaps otlier systems. A one-matrix model describes pure (i.e., not coupled 
to matter) surfaces and tlie bosonic string in zero dimensions. Its partition function is 



where (j) is an N x N hermitian matrix, and V is some non-gaussian potential; in the 
remainder of this work, it will be chosen to be 



Multi-matrix models describe random surfaces coupled to matter, or a string propagating 
in some target space. Their partition functions are of the form 



where the coupling matrix F, which depends on a matter coupling constant a, determines 
the matter model coupled to the surface, or the target space of the string. In this work, 
the free energy will always be defined in the spherical limit ^ oo. 

The matrix F may be thought of as the adjacency matrix of a labeled graph Q, which 
we shall call the target graph. One-matrix models have been solved exactly for a variety 
of potentials l^((/>) and for particular genera of the surface, as well as in the double scaling 
limit. Multi- matrix models have also been solved exactly, but only when the target graph 
^ is a tree, that is, when it has no cycles. This excludes a great many matter models (or 
target spaces), and in particular all those with central charge greater than one. 

A method to deal with arbitrary target graphs has recently been proposed 0: a 
low temperature expansion of the free energy, in powers of the matter coupling constant 
a. This reduces scalar multi-matrix model averages to contractions of tensor one-matrix 
model averages. Analogous to the low temperature expansion for ordinary lattice models, 
in n-th order one has some number of "blobs" which are open surfaces of uniform spin, 
connected by n edges which join blobs of unequal spin. 

In the present work, we consider regular target graphs, i.e., ones whose every ver- 
tex has the same degree or number of neighbors, A. We show, to all orders in the low 
temperature expansion, that in the limit A ^ oo the free energy of a multi-matrix model 
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becomes independent of its target graph and equal to some universal function, F^o, up 
to some simple transformations of a. Furthermore, we show that the only surfaces which 
contribute to the universal model are blobs joined into trees. Concretely, non-universal and 
non-tree terms in the free energy are suppressed by inverse powers of A. Some examples of 
models having regular target graphs, which we treat in detail, are the g-state Potts model, 
(i- dimensional lattice, and u copies of any regular model, where q, d, and u —>■ oo. Finally, 
we derive a simple, closed expression for Foo in any order. 

The d ^ oo limit of random surfaces has been discussed in the context of the 

continuum (i-dimensional model with Gaussian interaction for the embedding E Both of 
these works argue that that the worldsheet has a tree-like (or "branched polymer" ) struc- 
ture. Our results agree on this basic point. Beyond that, however, there are differences. 
Whereas we obtain entire open surfaces glued together into trees, according to [^j-Q the 
Gaussian model is dominated in this limit by worldsheet tubes connected into trees. There 
is no contradiction, because the models are different. Indeed, it is encouraging that some 
kind of trees emerge in both models. Furthermore, the "branched polymers" of appear 
at one of the critical regimes of our universal model (see last section): the critical point 
of the "matter" — it is to be remembered that the Gaussian model is always critical. The 
correspondence between our results and those of leads one to suspect that our results 
are valid non-perturbatively (beyond the low temperature expansion). 

This article is organized as follows. In the second section we summarize the old results 
on the low temperature expansion, and explain some new ones, in particular the explicit 
combinatorial rules for the diagrams in the partition function and in the free energy. In 
the third section we prove the assertions concerning universality and trees, and derive the 
exact form of F^o, which involves a sum over labeled trees. In the fourth section we go 
through some combinatorial manipulations to put each term in the sum into an algebraic 
form. In the final section, we discuss some of the implications of these results, including 
the critical behavior of the universal model, and suggest directions for further work. 

2. The low temperature expansion 

2.1. The target graph 

If the matrix model represents a bosonic string, the target graph Q is the discrete 
target space of the string, its edges the "metric." For instance, it might be reasonable 

^ To be contrasted with the Feynman propagators of the Cd model considered here 
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to take Q = Cd, the infinite (i- dimensional hypercubic lattice with periodic boundary 
conditions whose central charge is d. If the matrix model represents a random surface, the 
vertices of the target graph are the states of the matter, and the (weighted) edges are the 
interactions between the states.i For example, if the matter is a g-state Potts model, Q is 
/Cg, the complete graph on q vertices. The central charge is 0, ^, |, and 1 for g = 1, 2, 3, 
and 4, and undefined for q > 4 where KPZ scaling breaks down. So far, all edges of the 
target graph have the same weight, a. In either case, we shall use the word "spin" to 
designate the value of the vertex of Q assigned to each plaquette of the surface. 

In the case of a random surface with matter, one might also consider multiple "un- 
coupled" matter models [||. For instance, u species of, say. Potts models live on each 
elementary polygon of the surface, each spin interacting with nearest neighbor spins of the 
same species only. On a fixed surface this would be trivial: the partition function would 
just factorize. On a random surface, however, the different species effectively interact with 
each other because each interacts with the surface. The central charge of a multiple model 
is the sum of the individual central charges; the coupling matrix F is the direct product of 
the individual coupling matrices. As a result, the edges of the target graph of a multiple 
model never have uniform weights. For u Ising models, for example, the vertices of Q are 
those of a i^- dimensional hypercube, labeled 0, 1, . . . , 2^ — 1, but each vertex is connected 
to every other; any two vertices are connected with the weight a^, where 5 is the number 
of digits by which the binary representations of the labels of the two vertices differ. We 
shall call this target graph Q^. 

It should be noted that all target graphs mentioned above are regular. The degree 
A = d and — 1 for the lattice and Potts model, respectively. For multiple models, we 
define A to be the degree of the graph Q with all a^, a^, etc. edges removed; for u Ising 
models, therefore, Q is the ordinary (i-dimensional hypercube and A = u. In general, for 
u copies of a single model with degree 6, A = i/5. In this work, we shall treat only regular 
target graphs, although one would get similar results if this condition were to be somewhat 
relaxed. 

^ Strictly speaking, the "target graph" here should have adjacency matrix (1 — F)^^, the matter 
propagator. For the types of matter which we consider here — (multiple) Potts models — the two 
are equivalent up to a rescaling of the matrix fields and a transformation a — —f{a), where / is 
a monotonically increasing function in the relevant interval, and /(a) ~ a when a <C 1. 
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2.2. The partition function 

For a particular target graph Q with ^ vertices, the partition function of the matrix 
model is defined 

^^^^'"^ "\ / (2.1) 

The coupling matrix F depends on the matter coupling constant a (a = e~^, where /3 is the 
inverse of the matter temperature). We wish to expand the partition function in powers 
of a: 

oo 

^a(^,a) = l + E^«" (2.2) 



n! 

n=l 



For single models — where all edges of G are weighted or a — the n-th order coefficient 
Zn is a sum of possible terms of the form (Tr0^^0j^ • • • Tr where some of the 

indices may be equal. For a multiple model, in n-th order one has a sum of terms 
(Tr0j^0j^ • ■ ■TT(f)i^(f)j^) where m<n and 

m 

n ^Hi. = (2.3) 

k=l 

Each term {TrcpiCpj ■ ■ •) may be thought of as a subgraph H Q Q: the vertices of V, are 
a subset of those of Q, and the edges likewise, except that Ti, is allowed to have multiple 
edges. Each factor Tr 4>i4>j gives an edge in Ti between vertices i and j. Many subgraphs are 
isomorphic and therefore equal, for example (Tr(/)i(^2) and (Tr (/)2(/>3). We shall therefore 
take the subgraphs Ti to have unlabeled edges and vertices, because it is much more 
efficient; but in counting them, we shall remember that their edges are labeled. 

Consider, for instance, the term (Tr(/)i(/)2 (Ti (^24'^)'^ Tr030i), which might arise in 
fourth or higher order in the partition function. Its corresponding (unlabeled) graph H. is 




(2.4) 

The coefficient multiplying this term counts the number of ways in which the graph H.., 
with labeled edges, may be embedded in Q. In the qf-state Potts model {Q = JCq), H occurs 
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only in 24 with the coefficient ^q{q — 1)^; in u Ising models {Q = Qu), with the coefficient 
— 1) in Z5 and with other coefficients in higher orders; and H does not occur in 
the lattice Cd- 

The interpretation in terms of surfaces is as follows: each vertex of a subgraph Ti. is an 
open surface, a "blob." All spins on this blob are equal, or "frozen" at the same value, and 
therefore the surface can be described by a pure gravity matrix model (see below). Spins on 
two neighboring blobs are unequal. It should be noted that a blob contains contributions 
from both connected and disconnected surfaces. For instance, one of the graphsi which 
contribute to ( p^ ) is 




Each blob has many internal edges which connect equal spins. There are only four edges 
("bridges," dotted in the drawing) which connect blobs of unequal spin. We shall call the 
subgraphs Ti. "Z-skeletons" ; the meat on each vertex of a Z-skeleton is a blob, an open 
surface with uniform spin. 

As a concrete example, we give the low temperature expansions of the partition func- 
tions to third order for three different models. The Z-skeletons are represented graphically, 

^ We use the words "graph" and "surface" interchangeably for two objects that are dual to 
one another. 
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as above. First, the partition function of the qf-state Potts model, Q = KLq. 



Z = l + 



+ 



1 

+ 6 



+ 6 



+ 6 



+ 24 



+ 18 



+ 180 



9\» 



+ 18 



+ II2) -|(9-l)(6g'-35g2 + 77g-60) 



+ 0{a^) 



(2.5) 



Now a d-dimensional hypercubic lattice, Q = Ca with length L (which will drop out 
later in the calculations) and periodic boundary conditions {^ = dL'^): 



Z= l+l. .a 

£ *<=^ + 2{2d - 1)£ . 



1 

+ 2 



.+ {e-M + i)e 



a 



1 

+ 6 



£•^• + 2(2^- l)(2d-2)^^+6(2cZ- 1)^^^ 



+ 6{2d - l)e + 6{2d -l){£-6d + 2)1 \ + 3(£ - 4d + 1)£ l_ 

+ - 12de + 31 + 40^2 - 24d + 4)e 1= 

• 

+ 0{a^) 
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(2.6) 
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Finally, a multiple model: v Ising models, Q = Qv'- 



V . 



Z = 1 + -2^ 
2 
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3z/ 



. + 2v{y - 1) _^ + ^(^2^^ - 2Z. + 1) — + ^(z. - 1) . 
+ z/(z/ - l)(z/ - 2) + 3z/(z/ - 1)2 . . . . + 3z/(z/ - 1) 



{y - \){rv - 6z. + 4) 7^ + ^(z. - 1)(2V - 4z/ + 2) 



+ -{i^^'v^ - 122^2^2 + ^TTv + 40z^2 - 48z^ + 16) 
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+ 3z/2(^z/-l) 



+ - 1)(2'^ - 4) — + ^(^ - 1)(^ - 2) 



2'^a 



a 



(2.7) 



^.5. Calculating the blobs 

How can we calculate the contribution of each Z-skeleton to the partition function? 
This can be done by expressing multi-matrix averages of traces in terms of matrix com- 
ponents, so that they can be decomposed into contractions of one-matrix tensor averages. 
For example, our friend {Tr (j)i(j)2 (T^ <p2<p3)'^ Tr</'3<^i) can be decomposed into 



In the latter expression, the averages are with respect to a "pure gravity" one-matrix 
model: {O) = J Dfj,{(f))0{(j))/ J Dfi{(f)). The three tensors in this expressions are just the 
three blobs of the Z-skeleton A. Thus each blob with 6 adjoining bridges is a rank (5, 6) 
tensor, and each bridge contracts an upper and a lower index.@ 

The task is to calculate the one-matrix model vertex factor arising from a 5-valent 
blob: 

(2.9) 



It is necessary to distinguish between upper and lower indices to preserve the symmetry 
under unitary transformations. 



A convenient way to do this is to contract (|2.9| ) with A^^^ ■ ■ • A^^ , where A is some N x N 
hermitian matrix. Then the right-hand side of ( |2.9| ) becomes ((Tr A^)^), for which the 
generating functional is the external field integral 

Q{A,g) = J D^{ci^)e^^^^^ (2.10) 

P = log Q has been explicitly calculated in the spherical limit for an arbitrary matrix 
A by Kazakov and Kostov and by Gross and Newman . See for a summary of the 
result, as well as for explicit details of the vertex factor calculation. 

To understand the expansion of Q in powers of A, we put A ^ eA and expand P: 

P ^ Po - Pie + lrP2e^ - Ipse^ + ■ ■ ■ (2.11) 

Z D 

The first few terms are: 

pi = iV-^TrA {3g + 108g^ + lHQg^ + ■ ■ ■) 

= iV"^TrApi,i 
P2 = A^-^TrA^ {27g^ + 19Ug^ + ■ ■ ■) 

+ N-^{TrAf (9/ + 1296/ + ■ ■ ■) 
= N-^TrA^ P2,i + N-\TrAf p2,2 (2.12) 
Ps = A^-^TrA^ {6g + 540g^ + 58320£?^ + ■ ■ ■) 

+ A^-^TrA TrA^ (32^ + 69984£?^ H ) 

+ A^-3(TrA)3 {7776g^ + ■ ■ ■) 
= A^-^TrA^ ps,i + A-^TrA TrA^ ps 2 + A-^(TrA)^ ^3,3 

We observe that Pn,i/{n — 1)!, the coefficient of A Tr A" in p^ is just the connected n-point 
Green's function of the cubic one-matrix model @, which is as it should be. What, then, 
are the other terms which occur: Pn,2, etc.? The matrix A is attached to the external 
vertices, so one way to get, say, (TrA)^ is through a disconnected two-point function. But 
we know that P must be connected. The only other possibility is to have twisted graphs. 
For instance, a typical graph in p2,i is 



(2.13) 



while a typical graph in ^2,2 is 




(TrA)^ 



(2.14) 



Both of these graphs are "spherical," because by attaching tadpole endcaps to the external 
lines one obtains closed spherical graphs. But by tracing the indices in graph ( ^.141 ), one 
finds that the two A's never connect, giving the factor (TrA)^. The same holds for all 
Pn,T>i'- the surfaces which contribute to them can be transformed into ordinary connected 
Green's functions by untwisting the ends. 

The n-th order term in Q has the same contributions as its counterpart in P, the 
twisted and untwisted connected graphs. It also has, of course, twisted and untwisted 
disconnected graphs. 



2.4- The free energy 

The partition function looks rather complicated: it contains connected and discon- 
nected Z-skeletons, which contain connected and disconnected, twisted and untwisted 
blobs. If some of the blobs are disconnected, the whole graph may still be connected; 
but if the Z-skeleton is disconnected, the graph will be as well. Things are simplified 
immensely when one takes the logarithm to calculate the free energy 



Fg{g,a) 



1 



log 



/ia+^/2a2 + i/3a3 + .. 



(2.15) 



(we recall that ^ is the number of vertices ("volume") in the target graph Q). We can 
further simplify the calculations by expressing everything in terms of the coefficients Prn,T 
in the connected generating functional -P(A), rather than those in the actual disconnected 
blob-generating functional <5(A). 

When passing from the partition function to the free energy, one is of course left with 
only the connected graphs. Another major change occurs in the quality of the skeletons. 

Consider, for instance, one Ising model, whose target graph Q = Qi = • .. The Z- 

skeleton • — • — . does not occur in any order, because it is impossible to embed it in Q. The 
Z-skeleton however does occur in second order. Its bivalent blobs may be connected 
or disconnected, which gives rise to the following possibilities: 



^2.16) 
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The first of these is disconnected, so it gets canceled. The last is the skeleton The 
skeleton in the middle, however, is the graph • — . — ., which does not occur as a Z-skeleton. 
When we express everything in terms of connected blobs rather than the "natural blobs" 
which were first encountered, we find that many new skeletons are allowed in addition to 
the Z-skeletons, and that the combinatorial rules for finding the coefficients of the skeletons 
in the free energy are of course different as well. We shall call these new skeletons "F- 
skeletons." 

First, consider single models. As a concrete example, here is the free energy to third 
order of the g-state Potts model: 



F=^iq-l)Pla + 



1 



1 



-(g- 1) PiP2 + ^P2,l 



+ 



l(g - ifplPl +^{q- ifplPs + l(g - l)(g - 2)pl^ 



liq- 1)^2,1^3,1 + \iq- 1)W2,1P3,2 1)P3,1 



24' 



(2.17) 



and of the (i-dimensional lattice {D = 2d): 

F = ^Dpla + (^^D^PlP2 + \dpI,^ 

+ {^D^Plpl + ^D^PlP^ + ]^D'^piP2,iPz,i + \d'^PiP2,iPz,2 + ^^^^3,i) 0^ + C'(a^) 

(2.18) 

where pn = YliTPn,T- We would like to identify the above terms with F-skeletons. This is 
not hard, because a factor means that there are k n-valent vertices in the F-skeleton; 
r — 1 is the number of twists of the corresponding blob. This "degree sequence" uniquely 
identifies small graphs, so that we can make the following dictionary (omitting the twist 
index) : 



p\ 



PlP2 



PlP3 



p1 



P1P2P3 



2 2 
P1P2 



pI 



We are now ready to give the "Feynman rules" for the F-skeletons of single models. 
The F-skeletons Ti, are no longer subgraphs of the target graph ^, as we showed above. 
Instead, the vertices of Q are to be thought of as "colors" with which the vertices of V, are 
colored — the ordinary notion of spin. Two adjacent vertices of Ti. must have different 
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colors, and these two colors must be adjacent in Q. Thus the target graph Q provides 
both the possible colors for the F-skeletons, as well as the rules on how the colors may be 
combined. 

For a particular a graph H may appear as an F-skeleton provided that the latter may 
be colored according to the rules of the former. Consider some particular such unlabeled 
graph Ti. with n edges, and which therefore appears in the n-th order term f^. The 
coefficient of H in fn is a product of 

• the number of distinct ways to label the edges of Ti, which is nl divided by the order 
of the automorphism group of TC; 

• and a factor Cg{T-C), the number of ways to color each labeled version of H. 

Since /„ appears in the free energy divided by n!, and the free energy is divided by ^ 
(the volume of ^), the total factor of Ti. in the free energy is 

However, an F-skeleton 7i is represented in the free energy by more than one combination 
of the vertex factors Pn^r because the blobs may be twisted in different ways. For example, 
•=!>• — • may be represented by piP2,iP3,i or by piP2, 2^3,2, among others. Each of these com- 
binations may have a different coefficient, for which the expression ( |2.19| ) is an upper bound. 
The reason is topological: the whole graph must be spherical, while certain contractions of 
twisted blobs may yield non-spherical graphs. If one takes the untwisted n-valent vertex 
factors Pn^i for all vertices, one obtains the largest possible coefficient, which may still be 
less than the bound (|2T9|).i To illustrate for the F-skeleton — . in (|2T7| ) and ( pT8|) : 
the upper bound given by (|2.19| ) is ^{q ~ 1)^; this is saturated by PiP2,iP3,i] PiP2,iP3,2 is 
suppressed to ^{q — 1)^; the other four combinations are completely suppressed. 

The free energy rules are similar for multiple models, except that in order n there are 
F-skeletons with n and fewer edges, as there are edges in Q weighted a'^, a^, etc. For u 
Ising models {Q = Qu), the free energy to third order is 
1 



F = 2^Pia+ 



+ 
+ 



^iy^plp2 - \t^pI,1 + + 



^U^pIpI + ^I^^pIps + ^^^PlP2,lP3,l + ^^^PlP2,lP3,2 + ^^P3,l + 



(2.20) 



+ C(a*) 



This occurs for F-skeletons with cycles or multiple edges; see below. 
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For 1/ = 1 this can be checked against the exact solution for one Ising model with cubic 
vertices 0. One uses the same dictionary as before to translate from symbols to graphs. 
In each order n, the graphs with n edges follow the same combinatorial rules as single 
models. (To be more precise, the rules that would obtain for target graph Q, the graph 
Q with all a^, a^, etc. edges deleted. In the case G = Q^, G the ordinary z/-dimensional 
hypercube.) The graphs with fewer than n edges have modified coefficients, due to 

• a modified coloring factor CgiTi.n) (given a vertex of Q^, for example, there are p 
other vertices reachable by edges weighted a, (2) vertices reachable by a^-edges, and 
so on); 

• an "weight-entropy" factor which counts the number of ways that the weights of the 
edges can be chosen so that their weights satisfy (p73|). 

In a multiple model, once an F-skeleton appears in n-th order, it appears in all higher 
orders as well. The combinatorics of multiple models will be discussed in greater detail 
below. 



3. Large target graphs 



On examining the expansions (|2.17|) , (|2.18|) , and (|2.20| ) for the free energy of three 
different models, one notices some patterns: 

(a) the n-th order term is a degree n polynomial in A, the degree of ^ or ^; 

(b) the coefficient of A^c"^ in the free energy — which we shall call — contains contribu- 
tions only from those F-skeletons which are trees;! 

(c) for the single models (|2.17| ) and ( |2.18| ), the leading coefficients are, in fact, identical; 
the hn of the multiple model ( |2.2C1| ) are equal to the above if one ignores all F-skeletons 
with fewer than n edges; 

(d) the coefficients of the trees do not care about the twists of the blobs; in other words, 
the vertex factors pn,T in the tree coefficients always occur as the sums pn- 

These patterns are, in fact, not coincidental, and can easily be shown to hold for all 
regular target graphs and to all orders n in the low temperature expansion. Because of 
assertion (a) (see below), when A ^ 00 we can ignore all terms in order n but hn- Here 
we demonstrate the above assertions, and in the next section we derive simple expressions 
for h„. 



^ This is reminiscent of, but not identical to, the result of Erdos and Renyi [ p^ that finite 
random subgraphs of /d for n ^ 00 are almost always forests of trees. 
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The number of ways to label the edges of an F-skeleton Ti (which has v vertices) 
depends only on Ti and not on the target graph Q. For single models, the entire dependence 
of the coefficient of 7i on ^ is contained in the coloring factor CgiTi); for multiple models, 
the dependence on Q is in the modified coloring factor Cg (7i , n) . It is easy to see that 
CgiTi) has the upper bound 

Cg{H)<i/\'"-\ (3.1) 

which can be shown by picking an arbitrary vertex of Ti that can be colored ^ different ways, 
and thereafter traveling to every other vertex of 7Y, noticing that when one moves from a 
colored to a yet-to-be-colored vertex, one has at most A new colors to choose from. By 
induction, trees saturate the bound (|3.1|). The argument is analogous for multiple models. 
One has ^ colors for the first vertex of Ti, and for each subsequent vertex, traveling along 
an a'^-edge of one has at most ('^) < A'^ colors to choose from; using (|2.3| ), we obtain 
the same bound (|3.1| ). This demonstrates assertion (a). 

By "tree," we mean an F-skeleton which contains neither cycles nor multiple edges. 
Assertion (b) is now trivial, following from the upper bound ( |3.1| ) and the fact that a tree 
with n edges always has n -|- 1 vertices, while a connected non-tree graph with n edges 
always has n or fewer vertices. 

Having proved (b), we can ignore all graphs but trees when calculating h^- Since trees 
saturate the bound ( |3.1|) , their coloring factors are asymptotically (for single models) 

= Al^l+0(Al^l-i) (3.2) 

Because the coefficient of leading term is unity, the maximum contribution of a tree Ti to 
hni from ( p.l9| ), is 

^ (3.3) 



|Aut(7i:)| 

This maximum is, in fact, always saturated, for the following simple reason: one is free to 
twist any branch of a tree, without effecting the relative orientation of the other branches. 
In this way, one can untwist any twisted blob: all twists therefore contribute equally to the 
spherical average. We have demonstrated assertions (c) and (d) for single models, showing 
that the coefficient in of any tree F-skeleton is precisely (|3.2|). Assertion (d) is just as 
true for multiple models as it is for single models. Assertion (c) holds as well, provided 
that one restricts oneself to trees with n edges in order n, because when Ti has n edges, 
CgiTi.n) = CgiTi) and the weight-entropy factor is one. 
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As has already been mentioned, all four assertions are substantiated by the third-order 
expansions for the free energy for the Potts model ( p.l7|) , the lattice (|2.18|) , and multiple 
Ising models ( |2.2C1|) . The only F-skeletons which contribute to h^, for instance, are the trees 

• — • — • — • and ^/^- Their automorphism groups have order 2 and 6, respectively, which are 
the inverses of their coefficients for all three models. 



4. Counting trees 

We have shown that hn, the coefficient of A"^ in the n-th order term in the free energy, 
for any regular target graph describing a single model, is given by 

hn= y ^y^'^L, (4.1) 

^^JAut(T)| ^ ^ 

where Tn is the set of all unlabeled trees with n edges, and 

V{T)=p1^pI-... (4.2) 

where Vk is the number of /c-valent vertices in the tree T. By putting a b = aA, we 
expand F in inverse powers of A: F = Foo + 0(A~^), where Foo = Yl'^=i ^nb^. Given 
an unlabeled tree T with n edges and n + 1 vertices, there are (n + l)!/|Aut(T)| distinct 
ways of labeling its vertices. Therefore, -Foo can be expressed in a simpler way in terms of 
vertex-labeled trees: 



P'^-lEliEnT) (4.3) 

n=2 ■ re5„ 

where 5"^ is the set of all vertex-labeled trees with n vertices (of which there are n"~'^, 
by Cayley's theorem). Once again: F^o — the free energy of the "universal model" — is 
independent of the target graph, and the same for any single model. 

for multiple models can be expressed just as easily. First we derive it for u Ising 
models. A tree F-skeleton T with m edges first appears in order 6"^, and in every order 
thereafter. In m-th order, the coefficient of T is as for a single model; in higher order n 
it is modified, as explained above. Specifically, two neighboring points in T may have an 
edge in Q that is weighted a'^, a"^, etc. When coloring a new vertex in T along an a'^-edge, 
there are (^) colors to choose from, as opposed to u colors for single models. Since we 
are interested in the leading term in z^, we can approximate by u^/kl. The factors 
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u'^ combine to u"' in n-th order (see (|2.3|) ). But there are different combinations of edge 
weights which give z/"^: the sum is over all sets of m distinguishable positive integers whose 
sum is 71, where each term is the inverse of the product of the factorials of the integers. 
The reader can work out a few examples for himself, but it is not hard to see that the 
extra factor in n-th order, besides the usual l/|Aut(T)| term, is just the coefficient of 
in the Taylor expansion of (e^ — l)'", which is equal to 



k=0 ^ ^ 



(4.4) 



One obtains the same sum by transforming the coupling constant h in (^ 



1 



(4.5) 



Therefore, up to the transformation (|4.5| ) , F^o for multiple Ising models is equal to for 
single models. Note that now F^o is the coefficient of the leading power of v. The above ar- 
gument is easily generalizable to v copies of any single degree-A model; the transformation 
( ^75|) becomes 

b^b' = e^^- 1. (4.6) 



The sum over trees Xlres ^(^) been expressed in an algebraic form by Rior- 
dan[0. One defines the multi- variable Bell polynomials 



Yniyi, . . . ,2/n) = e" 



Jn 



x=0 



Riordan then shows that 



k=i 



Vk k 



(4.7) 



res„ 



(4.^ 



where — at the end of the calculation — one puts 6^ {n)k = n(n — l)---(n — /c + 1). For 
our purposes, this result can be re-expressed in a more convenient form by noticing that 
the Bell polynomial with the symbol 6 in every argument can be written 



Qn 

dx'' 
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[l + y{x) 



(4.9) 



x=0 



Therefore the sum over trees reduces to 



res„ 



d 



n-2 



dX 



n-2 



an 
ax 



A=0 



n(A, g)=p^x + ]-p2\^ + \p^\^ + 



(4.10) 



The vertex generating function n(A) is similar to the external field integral -P(A) 
(which generates the twisted vertex factors Pn,r, see ( p.lO|) ), except that it does not dis- 
criminate the vertex factors according to twist. Therefore it is given by the much simpler 
scalar version of -P(A): 



This may be calculated by the method of orthogonal polynomials, which gives 



(4.11) 



n(A,^7) 







dx (1 — x) log ^ ^ 

X 



1 1 Kg) 



(4.12) 



where the function r satisfies the equations 

—Qgr = —A + s + ?>gs^ 

-6gX = (1 + Qgs){-X + s + 3gs'^) 



(4.13) 



In order A"^ and higher, the leading term in g in 



and h{g) = dxexp ^/NXx — V{x 
n(A, g) comes from the second term in ( [1.12|) , and in order A"^"'"^ is equal to 

(2n)! 



n!(n + 2) 



■(3^)' 



(4.14) 



5. Discussion 

We have shown that for regular target graphs of degree A the leading term of the free 
energy in inverse powers of A, Foo, is identical (for single models) for all target graphs and 
contains contributions only from trees. In particular, after rescaling the matter coupling 
constant a b = a A, the leading term is given simply by equation ( [1.3|) . This, after some 
combinatorial manipulations and simple observations concerning trees, can be shown to 
be equal to 



d 



n-2 



n=2 



an 

aA 



(5.1) 



A=0 
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where n(A, g) is the one-matrix model average ( [4.11|) . An almost identical result holds for 
multiple models. For u copies of a model with target graph Q with degree A, the leading 
term in the expansion of the free energy in inverse powers of u is identical to (|5.1| ) after 
the transformations a ^ b = au and b ^ b' = e^^ — 1. These results have been shown to 
hold in every order of the low temperature expansion. 

If A becomes very large — and therefore the target graph becomes very large — we can 
ignore all but the leading term Foo- For the lattice and multiple 2, 3, or 4-state Potts 
models, this limit corresponds to central charge c — > oo. The random surface picture is 
that of a series of "baby universes" (blobs) connected to each other as a tree (F-skeleton). 
The spins on each baby universe are frozen equal, while the spins on neighboring baby 
universes are unequal. 

For ordinary multi-matrix models (on small target graphs), the cosmological constant 
g is coupled to the area of the surface, while the matter coupling constant a is coupled 
to the energy of the matter on the surface. The typical F-skeletons occurring in these 
models are highly multiply connected graphs. There are three critical regimes: adjusting 
g makes the area blow up; adjusting a orders the spins; and adjusting both yields a critical 
matter model on a continuous surface. For large target graphs, when we are left with the 
universal model which has only trees, the picture is different. The cosmological constant 
is still coupled to the area, but here this is the area of each baby universe. The matter 
coupling constant, however, is now coupled to the size (the number of edges, and therefore 
the number of vertices) of the tree. Adjusting g makes the areas of the baby universes 
blow up. Adjusting a makes the size of the tree — the number of baby universes — blow up. 
Finally, adjusting both yields an infinite tree whose vertices are continuous surfaces, the 
baby universes. 

So one still gets a phase transition by adjusting the matter coupling constant, only it is 
a tree-growing — rather than a matter-ordering — transition. The matter coupling constant 
becomes a "second-quantized" cosmological constant. The critical exponent, a, of the free 
energy is likely to be ^. The reason is that if one puts the vertex factors = 1, the sum 
over trees just reduces to the Cayley value n"^"^ in n-th order; the coefficient of 6"^ grows 
as 71 """^/n! ~ n~^/^e"', which gives a = |. The n~^l'^ behavior is robust: it seems to be 
independent of the values of the vertex factors p„; we therefore expect that a = \. It 
would be interesting to investigate the critical behavior of Foo systematically. 

Interestingly, \ is the value of the exponent found in and [Q, but for the "first- 
quantized" exponent 7str rather than the "second-quantized" exponent a. For a gaussian 
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model, of course, there is only one coupling constant and — generically — only one critical 
regime. Comparing our results with those of 0, we see that for the gaussian models 
one can blow up the trees, but not the baby universes. Perhaps this deviation from the 
universal model indicates a subtlety in the target space continuum limit. 

Finally, we would like to point out the resemblance between our results and mean 
field theory. In ordinary spin models, the free energy becomes universal (i.e., independent 
of the type of matter) when the dimension of the worldsheet gets large. Here, the free 
energy for a two-dimensional worldsheet is universal when the dimension (or number of 
nearest neighbors) of the target graph is infinite. In ordinary spin models, the transition 
to universal, mean-field behavior occurs abruptly at a finite dimension. Could this also be 
the case for random surfaces? 
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